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Abstract 

A simple formula for computing the generalized Hubbell radiation rectangular source integral 
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is introduced. Tables are given to compare the numerical values derived from our approximation 
formula with those given earlier in the literature. 
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1 Introduction 

In their pioneering work, Hubbell et al. (1960) obtained a series expansion for the calculation of 
radiation field generated by a plane isotropic rectangular source (plaque), in which the leading term 
is the integral 

. . (J f ( a \ dx , . 

Iia.b) = — / arctan - , (1) 

here a is the uniform surface source strength per unit source area. In equation (TIJ the quantities 
a = w/h and b = l/h are defined in the range < a < b < 00, where h is the height over the 
a corner of a plaque of length I and width w. For the important applications of this integral in 
many problems in radiation field, different methods were introduced to obtain numerical values of 
detector response to plaque source 
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For some of these methods we refer to the work of Glasser, 1984; Kalla et al., 1987; Ghose et al., 
1988; Galue et al., 1988; Gotze, 1995; Kalla, 1993; Timus, 1993; Michieli and Maximovic, 1996; 
Kalla and Khajah (1997, 2000); Prabha, 2001; Stalker, 2001; Guseinov et al., 2004; Ezure, 2005; 
Guseinov et al, 2005; Prabha, 2007. 

Although h(a,b) is not expressible in simple closed form, Glasser (1984) has evaluated it in 
terms of Appell's hypergeometric function F 2 (see for example Slater (1966), Ch. 8, for a study of 
Appell functions F q , q = 1,2,3,4). Indeed, elementary differentiation of ([2]), with respect to a, we 
have 

/o l + x 2 + y 2 { ' 

and straightforward substitutions x = a^/u and y = byjv allow as to write ^ as the double- 
integral representation of Appell's hypergeometric function F 2 (Slater (1966), Ch. 8, formula 
8.2.3), consequently 
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Various generalizations of equation (TjQ) have been given in the literature (see for example, Kalla et 
al., 1987; Galue et al, 1988; Saigo and Srivastava, 1990; Galue, 1991; Kalla, 1993; Galue et al, 1994, 
Kalla et al. 2002, Oner 2007). More specifically, Kalla et al. (1987) introduced a generalization 
defined by the integral 
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where 7 > j3 > 0; a, b,p > 0; — 1 < A < 2a — 1; and 2-^1(0, P) 7! x) is Gauss hypergeometric function 
(Slater, 1966, Ch. 1). We notice that 
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by virtue of the identity (Slater, 1966; formula 1.5.11) 
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By selecting suitable values for the parameters a, f3 and 7, equation ([5|) can be reduced to different 
integrals with potential applications in radiation-field problems of specific configurations of source, 
barrier and detector (Kalla, 1993). Such results are also useful in illumination and heat-exchange 
engineering Boast, 1942; Fano et al., 1959; Hubbell, 1960. Using a simple transformation, x = by/u, 
equation © can be written as 
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which is easily compared with the single-integral representation of the Appell hypergeometric func- 
tion F2 (Opps et al., 2005, formula (2.6)) to yields 
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Recently, Opps et al. (2009) establish a number of new recursion formulas for the Appell hypergeo- 
metric functions F 2 wherein some applications to the evaluation of some generalized radiation field 
integrals were discussed. The purpose of the present work is to continue our investigation of finding 
closed form and approximation formulas for effectively computing the radiation field integrals such 
as equations (PQ) and ([5]) . In the next section, we develop a new approximation formula to evaluate 
precisely and to any desire degree of accuracy the generalized Hubbell radiation rectangular source 
integral ([9]). In section 3, numerical results and comparisons with previously reported values are 
presented. 

2 The computation of F 2 (a- (3, A±l ; 7> A+3 ; _£) 

May be one of the most important cases regarding the computations of the radiation field integrals 
that capture the interest of many researchers is evaluating the integral Eq. (fTJ) effectively and 
precisely. Some researchers were able to evaluate h(a, b) using rapidly convergent series (see the 
original work of Hubbell et al, 1960), further Gabutti et al. (1991) investigated h(a,b) in terms 
of its series expansions while numerical computations of this integral have been carried out by 
Hanak and Cechak (1978), Gotze (1995) developed an effective method for computing the Hubbell 
radiation rectangular source integral, Kalla and Khajah (1997) (see also Kalla and Khajah (2000)) 
used Tau Method to approximate H(a, b), Stalker (2001) used new convergent series for evaluating 
h(a, b) for large a and b and more recently, Ezure (2005) used Haselogrove method, Guseinov et al. 
(2004) used binomial expansion (see also Guseinov et al. (2005)), Prabha (2006) expressed again 
h(a,b) using some recurrence relations (see also Prabha (2007)). For a survey of various methods 
in computing the Hubbell rectangular source integral Eq.(fT|) and its generalization, we refer to the 
work of Kalla et al. 2002. 

In this section, we given a new approximation equation that can be used to compute H[ a i ^'^'^'] 
to any degree of precision and, byproduct, we can, therefore, evaluate the Hubbell radiation rectan- 
gular source integral ([T]). Our approximation expression based on the following recurrence formula 
for F 2 (see Opps et al. (2009) for detailed proof.). 

Theorem 1: For \x\ + \y\ < 1; n > 0; a, «i, a 2 6 C; /3 2 G C\Zq , the Appell hypergeometric 
function F2 satisfies the following identity 

F 2 (o-,ai,a<2-n;/3i,(3 2 ;x,y) = F 2 (a;a 1 ,a 2 ; Pi, fa;x,y) 

n 

~ ^y)i J 2(ff + l;Q!i,a2-*; + l;/3i,/32 + l;a;,y). □ (10) 

Writing F 2 (a;/3,^±i; 7,^; -£,-£) as F 2 (a;/3,^- 1; 7, -f, -f) and apply the recur- 
rence relation (I10D. we obtain 
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we may now write Eq. (|ll|) as 
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Further, we may now regard F 2 (a+1 ; /3, ^ ; 7 , ^ ; - f , - f ) as F 2 (a+1; /3, ^ - 1; 7 , ^ ; - ^ , - f ) 
and apply the recurrence relation (jlOp again to obtain 
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After similar n — 2 steps, we arrive at 

2 2 p p ^2±Aj^ p p p + b z 
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(15) 

where (ce)k denotes the Pochhammer symbol defined, in terms of Gamma functions, by 

__ r(q + k) _ f 1 if (fc = 0; a € C\{0}) , 

(ajfe ~ r(a) ~ \ a(a + l)(a + 2)...(a + k-l) if (jfc e N; a e C), 

here N being the set of positive integers. 

In principle, the computation of F 2 {a + n + 1; (3, A+2 2 n+3 ; 7 , A+2 2 "+ 5 ; — — y) follows the same 
technique and consequently, for large n 

F 2 a;/3,— -; 7 , __; - — « £ -A-^- - 1 + - *F X (a + k, /3; 7 ; — ^ • 

2 2 p p ^(3±A) fc V;p/ \ pj p + b 2 

(16) 

From which we now have for large n, 

since the two-terms asymptotic expansion of the Appell hypergeometric function F 2 developed by 
Lopez and Pagola (2008) indicate that for large n, the Appell hypergeometric function 

n A + 2n + 3 A + 2n + 5 a 2 6 2 . 

F 2 (a + n+l;/3, ; 7 , ; , ) 

2 2 p p 

on the right-hand side of Eq . (|15j) approach zero. 
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3 Numerical results and discussion 



In order to test our approximation formula (|17p and to show that it indeed simplify much of the 
numerical complexity involved in calculating the radiation field integrals such as ([1]) and ©, we 
compare, first, our approximation formula against the exact equation obtained earlier in computing 
a, b, 1, 



H 



l l 

2> 2> 



1 



H 



(Opps et al. (2009), equation (83)), namely, for a = 1, 
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In Table 1, we reported our computation using equations ([IB]) and (fTUj) for different values of a, 
6 and p, the calculations were performed using MATHEMATICA software version 7. It should 
be clear that any discrepancies may appear are due to the numerical accuracy used in computing 
the Gauss hyper geometric functions. In Table 2, we compared our calculated values with those 



Table 1: Comparison of the values of H [ a { f' ] for some values of a, b and p calculated using 

2 ' 2 ' 



equations 


CHI) 


and CED 






a 


b 


P 


Eq. (pSD 


Eq. CE]) 


0.1 


0.2 


0.5 


0.000 219 698 305 361 162 27 


0.000 219 698 305 361 161 92 


0.1 


0.5 


0.5 


0.001 259 518 891 975 573 3 


0.001 259 518 891 975 572 7 


0.2 


0.2 


2.0 


0.000 222 868 191 957 851 78 


0.000 222 868 191 957 853 9 


0.2 


1.0 


2.0 


0.005 038 075 567 902 293 


0.005 038 075 567 902 291 


0.5 


0.5 


0.5 


0.005 794 884 270 704 959 5 


0.005 794 884 270 704 952 5 


0.5 


1.0 


2.5 


0.011 293 885 774 813 334 5 


0.011 293 885 774 813 332 



obtained by Guseinov and Memedov (2005) and with those obtained by Galue et al. (1994). It 
should be clear that equation (|17p can be used for arbitrary values of the parameters a, b, p and it 
is not restricted to any particular range of parameter values. 

In Table 3, we report our numerical computation of the Hubbell rectangular source integral ([TJ 
for a = p = 1 as well for some other values of p using equation (I17p , or simply 

H [ ? H'° 1 « ^ ± " (») 1 (l + + ,, i; | ,20) 

- x '2>2 J 4 ^p^(l)fcVp/ V pj 2 2 P + b 2 

In the same table, we also compared our results with the earlier numerical values obtained of 
Guseinov and Mademov (2005) and Galue et al (1994). It is worth noting an important feature of 
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Table 2: The values of H[ a {'i'}] integrals for a = 1 and some values of a, b and p obtained from 

2 ' 2 ' 

Eq. (H7D, Guesinov and Mamedov (2005) and Galue et al. (1994). 



a 


b 


p 


Eq. (HID 


Guseinov &; Mamedov 


Galue et al. 










(2005) 


(1994) 


0.1 


0.2 


0.5 


0.000 219 698 305 361 161 92 


0.000 219 698 305 352 979 


0.000 219 698 31 


0.1 


0.5 


0.5 


0.001 259 518 891 975 572 7 


0.001 259 518 892 096 95 


0.001 259 518 9 


0.2 


0.2 


2.0 


0.000 222 868 191 957 853 9 


0.000 222 868 191 569 00 


0.000 222 868 19 


0.2 


1.0 


2.0 


0.005 038 075 567 902 291 


0.005 038 075 568 393 22 


0.005 038 075 6 


0.5 


0.5 


0.5 


0.005 794 884 270 704 952 5 


0.005 794 884 271 955 31 


0.005 794 884 3 


0.5 


1.0 


2.5 


0.011 293 885 774 813 332 


0.011 293 885 774 293 


0.011 293 886 



Table 3: The comparative values of H[^'^ 5 ] integral from Eq. (|17j) . Guseinov and Mamedov 
(2005) and Galue et al (1994) 



a 


b 


p 


Eq. ([17]) 


Guseinov & Mamedov 


Galue et al. 










(2005) 


(1994) 


0.1 


0.1 


0.5 


0.001 570 716 369 171 686 


0.001 570 716 369 157 32 


0.001 570 716 4 


0.1 


0.5 


1.0 


0.003 678 199 808 681 331 


0.003 678 199 808 778 47 


0.003 678 199 8 


0.2 


0.5 


2.5 


0.003 067 148 756 523 266 


0.003 067 148 756 507 53 


0.003 067 148 8 


0.2 


0.8 


2.0 


0.005 758 600 701 534 181 


0.005 758 600 701 357 22 


0.005 758 600 7 


0.5 


0.5 


1.0 


0.017 188 506 077 049 23 


0.017 188 506 077 717 6 


0.017 188 506 


0.6 


0.5 


1.0 


0.020 067 469 440 496 68 


0.020 067 469 441 718 5 


0.020 067 469 


0.8 


0.6 


2.8 


0.012 248 693 964 171 94 


0.012 248 693 963 979 3 




1.0 


0.8 


4.2 


0.013 484 796 561 457 52 


0.013 484 796 561 864 6 




0.5 


2.0 


5.4 


0.012 012 547 384 013 146 


0.012 012 547 385 362 9 




0.8 


2.6 


7.5 


0.017 255 112 588 899 273 


0.017 255 128 853 560 6 





the approximation formula ()17p is that it is self-adjusting; that is, if 

\H n+1 -H n \£e (21) 

where e is the desire accuracy then n should be increased to reach the required accuracy. Here, H n 
refer to the right-hand side of equation (|17p . 
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